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DISSIPATIVE DYNAMICS IN SEMICONDUCTORS AT LOW 

TEMPERATURE 

GEORGE ANDROULAKIS, JEAN BELLISSARD, AND CHRISTIAN SADEL 

Abstract. A mathematical model is introduced which describes the dissipation of elec- 
trons in lightly doped semi-conductors. The dissipation operator is proved to be densely 
defined and positive and to generate a Markov semigroup of operators. The spectrum 
of the dissipation operator is studied and it is shown that zero is a simple eigenvalue, 
which makes the equilibrium state unique. Also it is shown that there is a gap between 
zero and the rest of its spectrum which makes the return to equilibrium exponentially 
fast in time. 



I~~5 The sciences do not try to explain, they hardly even try to interpret, they mainly make models. By a model 

\^ is meant a mathematical construct which, with the addition of certain verbal interpretations, describes 

^_^ observed phenomena. The justification of such a mathematical construct is solely and precisely that it is 

1-^ expected to work. (J. von Neumann '51') 

^ 1. Introduction 

, Q, This article is dedicated to the construction and the fundamental properties of a model of 

dissipative transport, describing the electron or hole transport in semiconductors at very 

low temperature. By "very low" it is meant that the temperature is low enough so as 

00 to confine the charge carriers to the impurity band. Without dissipation, the transport 

^ is coherent and is likely to be described by an Anderson model, namely a Schrodinger 

operator on a discrete lattice with a random potential. In lightly doped semiconductors 

this model has to be considered in the strong localization regime as will be explained 

O ill Section [2j The main source of dissipation in all solids, including semiconductors, is 

^~~' coming from the electron-phonon interaction. Namely, the Coulomb interaction between 

electrons and nuclei leads to the slow nuclei motion when an electron is passing by. The 

.!_I harmonic interaction between nuclei leads, in turn, to the production of acoustic waves, 

/\ that are quantized, at least if produced in small quantities. These quanta of acoustic 

j^ waves are called phonons. These waves propagate in the crystal at the speed of sound 

and can kick another electron eventually, leading to loss of information about the electron 

quantum state. There are other sources of dissipation for the electron motion like the 

direct Coulomb interaction between them, spin coupling or photon emission. However, 
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it has been shown [22] that the first process is quantitatively more important than all 
others. Still, even the dominant mechanism for dissipation can be considered as weak in 
most cases, in particular for the problem considered in the present paper. Therefore, it is 
legitimate to approximate the system by using the so-called Markov approximation. As 
a consequence, the dissipative dynamic can be described through a Markov semigroup 
whose generator is a Lindblad operator [28], also called nowadays Lindbladian. Like 
in Atomic Physics \i2i, this Lindbladian could be computed from their second order 
perturbation theory, called the Fermi golden rule. It will not be the method used here. 
The purpose of the present work will be to construct a Lindblad operator that describes 
phenomenologically the dominant sources of dissipation, and to investigate its spectral 



properties. The precise description of the model is given in Subsection 2.5 The spectral 
properties are summarized in Section |5} The main result is the Theorem 5.5, showing 
that (i) such an operator and the dynamical semigroup it generates are well defined even 
if the system is out of equilibrium and that (ii) it forces the return to equilibrium if there 
is no gradient of chemical potential or of temperature in the system. However, the main 
new contribution of the present work lies in the mathematical framework as explained 
below. But since the explicit model is strongly dependent of the physical regime that it 
intends to describe, it will be necessary to describe the physics in detail in order to make 
sure that the model is realistic (see Section [5]). 

The problem investigated here is motivated by the mechanism called variable range hop- 
ping conductivity [39]. It has dominated the study of semiconductors for almost two 
decades since the work of Miller and Abrahams in 1960 ^SlJ and the seminal contribution 
of Mott [32] predicting the behavior of the conductivity as a function of the temperature. 
It has been suggested [SB] that this regime is also the dominant contribution in the Quan- 
tum Hall Effect [5] that explains in particular the amazing accuracy of the experiment. 
For indeed, the relative error with which the Hall conductivity can be measured in a QHE 
experiment is of the order of 10"^" due to the smallness of the direct conductivity on 
the plateau of conductivity [34J . The theory of the integer QHE has been made rigorous 
through the use of the formalism of Non-commutative Geometry in [5] . In this latter work, 
the problem of dissipative transport was investigated within the so-called relaxation time 
approximation (RTA). The RTA reduces the dissipative mechanism to only one time scale 
and allows to consider the charge carriers as independent particles. As shown in [S], when 
applied to the QHE, this approximation leads to a relative error of 10~^ with the best 
data, namely six orders of magnitude larger than what is actually observed !! The rea- 
son is that the charge carriers conductivity is suppressed by the variable range hopping, 
as was shown by Mott [32]. However, such a mechanism involves an infinite number of 
time scales. It was advocated in [6] that such a mechanism can be represented through a 
Markov semigroup, the generator of which is a Lindblad operator. Unfortunately previous 
attempts to implement this idea have provided mixed results. The main reason is that 
the charge carriers, electrons or holes, are Fermions and any approximation leading to 
consider these particles in a semi classical regime fails to include the statistical correla- 
tion induced by the electron-phonon interactions. In order to successfully represent this 
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mechanism, a mathematical model must take second quantization into account. The main 
new contribution of the present paper is precisely to work with a many-body formalism. 
Since this approach is technically very demanding, the model will be simplified. The main 
simplification consists in replacing the Anderson model, describing the coherent part of 
the motion, by a purely potential contribution, neglecting the kinetic part, which, in real 
semiconductors is indeed extremely small. This kinetic part will be reintroduced in the 
dissipative mechanism through a contribution of the tunneling effect between impurities. 

Even with so many simplifications, the formalism is heavy and will occupy most of this 
paper. This is because the random character of the distribution of impurities breaks 
the translation invariance. Since the early eighties, thanks to using the ideas of Non- 
commutative Geometry [15], the formalism required to replace the translation group is 
known (for instance, see [6] and references therein): a groupoid replaces the group of 
translations. The inclusion of the many-body formalism in this framework was developed 
in the PhD Thesis of Dominique Spehner [H]. This leads to replace the observable 
algebra by a continuous field of C*-algebras over the space describing the disorder. The 
notion of continuous field of Banach spaces was introduced by Tomyama [171 HHl 09] in 
the context of the spectral theory for C*-algebras and later developed by Dixmier [T7] . 
While the concept is easy to understand, it is technically demanding. Then, the coherent 
dynamics can be defined as a continuous field of dynamics, leading to a continuous field of 
KMS-states describing the equilibrium of the electron gas in the solid. In much the same 
way, the dissipative dynamics is defined by a continuous field of Markov semigroups. In 
the present work, various existence results for the dynamics are proposed. One, mainly 



the Theorem 4.4, is based on the estimates used by Bratteli and Robinson [TU] to prove 



the existence of the dynamics in the many-body theory. The other one, mainly the 



Theorem |4.9[ uses the continuous field of Hilbert spaces generated by the field of KMS- 
states, through the GNS-construction, and sees the Lindbladian as a non-commutative 
analog of a Dirichlet form. Dirichlet forms were introduced by Beurling and Deny [8] and 
the characterization was completed by Fukushima [20]. The noncommutative Dirichlet 
forms were defined by Albeverio and Hoegh-Kr0hn [2] and they were characterized in 
full generality by Cipriani [11]. The definition requires some notion of positivity in the 
Hilbert space of the GNS-representation. Such a positivity is provided by a cone in the 
Hilbert space, that was identified and characterized in the early seventies by Araki [1] 
and Connes [T3|. In each fiber of the field of Hilbert spaces, it is the cone generated by 



positive elements of the corresponding fiber of the field of observable algebras. For the 



sake of the reader, this will be explained in Section 3^ 

Because the formalism required here is so heavy, it seems wiser to restrict the present 
paper to the description of the dynamics and the return to equilibrium. However, the real 
goal is to show that this model is liable to account for the Mott prediction concerning the 
low temperature behavior of the conductivity. An important result was obtained by [H] : 
by looking at the variable range hopping at very large length scale, as a random walk in a 
random environment the authors could prove that the Mott prediction was a lower bound 
to the diffusion constant. However, to get an upper bound is highly non trivial. In order to 
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do SO, using the present model, it will be necessary to face two challenges. The first one is 
the definition of the local currents. As it turns out, this is not a trivial problem because of 
possible divergent effects in the infinite volume limit. It requires insight about the physical 
nature of currents and of the dissipation mechanisms. With the proper concept, though 
it is possible to prove rigorously the validity of the Kubo formula [I] whenever the charge 
carriers can be considered as a continuous fluid. The other challenge is the discontinuous 
nature of the variable range hopping mechanism, forcing the charge carriers to hop at 
distances ten times larger than the average distance between impurities. In particular, 
the continuous fluid picture breaks down ! This is why Miller and Abrahams |31] proposed 
to see the charge carriers as electric currents in a random network. In the early seventies 
percolation theory was successfully introduced into this picture [3|, ESI EH] to justify the 
prediction of Mott. However, several approximations, justified by the physical situation, 
require additional work in order to make this argument rigorous within the scope of the 
present model. It will be hopefully the subject of a forthcoming pubhcation. 

Acknowledgments: This work benefited from the NSF grants DMS-0600956 and 
DMS-0901514. Part of this work was done in Bielefeld with the support of the SFB 701 
"Spectral Structures and Topological Methods in Mathematics" during the Summers 2009 
and 2010. G.A. and CS. thanks the School of Mathematics at the Georgia Institute of 
Technology for support during the Spring 2009. 



2. Physics of lightly doped semi-conductors 

2.1. Orders of Magnitude. The content of this section can be found in several textbook, 
in particular the one by Shklovskii and Efros |59| . 

The two types of semi-conductors that are the most used and studied today are silicon and 
GaAs, due to their importance in modern electronics, telecommunication and computer 
hardware. Silicon is currently obtained, in the industrial process, in a form of cylindrical 
ingots of about 2m in length and 25 cm in diameter. The crystal is perfect with less 
than 10~^° impurity or defect per atom. Because silicon has 4 valence electrons, the 
atomic orbitals hybridate in the sps form, leading to a diamond crystal, where tetrahedra 
alternating in a staggered way. Ga belongs to the column III of the periodic table, namely 
it has 3 valence electrons. As has 5 valence electrons and thus, it belongs to the column 
V. Mixed together in equal quantity, Ga and As exchange one electron to produce pairs 
of tetrahedra, leading also to a diamond lattice. If the electron-electron interactions are 
ignored, the band spectrum is similar to the one of graphene, namely two bands touch on 
a conical point exactly at the Fermi level, (the Fermi level is the maximum energy of an 
electron at zero temperature). 

However, the Coulomb interaction between the valence electrons belonging to the same 
atom, is strong enough to force the opening of a gap at the Fermi level, which will be called 
the main gap. For Si, this gap is about 1.1 e\^, which is enormous. Correspondingly, the 
temperature necessary to allow a large number of electrons to jump from the valence band 
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Figure 1. (i) Left: energy ranges (ii) Right: 1 & 2 hydrogen atom levels 



to the conduction band would be 1.2 x 10^ K . In other words, the clean silicon is a perfect 
insulator. For GaAs the gap is 1.43 e^ corresponding to a temperature of 1.6 x 10^ -ft'. 

For the purpose of electronic applications, the doping is about 10~^ impurity per atomic 
site. Such a crystal is called lightly doped. For ra-type doping, the impurities are atoms 
with one more electron in the valence band than the crystal atoms (such as P,As,Sb). 
Then, an impurity band is created in the gap at a distance of about lOmeV from the 
conduction band (see Fig. [l]). Hence, the gap between the impurity band and the conduc- 
tion band is 100 times smaller than the main gap. In particular, at room temperature, 
most electrons of the impurity band jump into the conduction band, which explains why 
semi-conductors are actually conductors. However, the conduction electron density is 
controlled by the impurity concentration, namely it is much smaller than in metals. At 
very low temperatures, say around 1 K, this gap is too large to allow electrons to jump, 
so that electrons are confined in the impurity band (see Fig. [I]). Similarly, for a p-type 
doping, the impurities are acceptors instead, namely they have one electron less on their 
valence band than the crystal atoms (such as Al,Ga,In). Then, a hole-impurity band 
is created near the valence band at approximately the same distance as before. Similar 
conclusions arise after proceeding to the hole-particle symmetry: (i) the charge carriers 
are holes instead of electrons, (ii) the origin of the energy is the top of the valence band 
instead of the bottom of the conduction band and (iii) the sign of energy is reversed. 

For an impurity concentration c ~ 10~^ per atom, the average distance between impur ities 
is given by Vav = c~^^^ ~ 1000 atomic spacings. As will be seen in the next Section 2.2 
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an isolated impurity behaves like an Hydrogen atom with a re-normalized mass and a re- 
normalized Coulomb coupling constant. In these materials, the Bohr radius of an impurity 
is ttB — lOOA. Hence the average distance between impurities corresponds approximately 
to lOaB- This gives a band width of approximately ly ~ 1 meV, 10 times smaller than 
the gap between the impurity band and the conduction (or valence) band (see Fig. [I]). 

In most cases a semi-conductor has the two types of impurities simultaneously. In such a 
case it is called compensate. If n-type impurities dominate, then all acceptors will get an 
extra electron since their energy is lower, creating negative ions, and some proportion of 
the donors will become ionized creating positive charges in the crystal. The Fermi level, 
namely the maximum energy of an electron at zero temperature, will be located within the 
donor impurity band at a position depending upon the relative concentration of donor and 
acceptor. If the acceptors dominate, some of them will acquire all the electrons coming 
form the donors and will become negatively charged, while all the donors will be ionized. 
The Fermi level will then be inside the acceptor impurity band and the conduction will 
be due to holes instead. In both cases however, the positions of the ions are random, 
creating a random Coulomb field within the crystal which influence the energy level of 
the impurity electrons. 

2.2. One Impurity &: Hydrogen Atom. If one donor impurity is inserted in a perfectly 
periodic crystal. Slater [40j showed in 1949 using a theorem by Wannier ^\, that the 
extra electron behaves like in an hydrogen atom, with re-normalized parameters for its 
mass and the dielectric constant, provided the origin of energies is taken at the bottom 
of the conduction band. If the impurity is an acceptor instead, a similar argument can 
be used for holes instead of electrons, using electron-hole symmetry. If the bottom of 
the conduction band is an isotropic minimum, the corresponding effective Hamiltonian 
describing the electron near the impurity is given by 

Hh = -2EA^ + Y-A , r = —,aB=^,Eo = —^. 2.1 

\ 2 \r\J as rne"^ Ih'^K'^ 

Here, e is the electron charge, m is the effective mass of the electron in the crystal and k 
is the re-normalized Coulomb constant in the medium. The ground state is given by 

0o(x) = ^^, Z,= [ixaiy'' (2.2) 

so that a^, called the Bohr radius, gives the length scale beyond which the wave function 



becomes negligible. The corresponding eigenvalue is —Eq given in eq. (2.1). The other 
energy levels are given by E^ = —Eq/ti^ for n = 1, 2, 3, ■ ■ ■ . In particular the gap between 
the ground-state and the first excited state is "iEo/A (see Fig. [T|. 

As a result, in the donor case, the lower energy state of the impurity is located in the 
gap of the crystal, at a distance Eq from the conduction band. Applied to impurities in 
a crystal, the Wannier equivalence theorem provides effective values for m and k, which 
are usually different from the values for the electron in the vacuum. For silicon, this Bohr 
radius is about lOOA and Eq ^ 10 meV [39] . 
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2.3. The Anderson model. In a semi-conductor, there is a density of impurities. Even 
if this density is as smaU as 10~^, as it is for hghtly doped media, the effect on the charge 
carrier dynamics is not negligible. The first important effect occurs if the semi-conductor 
is compensated. For indeed since the density of donors and acceptors are not equal, 
there is a nonzero density of ions in the crystal. The position of these ions is random. 
Consequently, they create a static electric field within the crystal, which is also random. 
Since the charge carriers are confined on the impurity site, they cannot move to screen 
this field. As a consequence, the effective potential seen by each charge carrier at an 
impurity site, is itself shifted by a random term. Since the distance between ions is very 
large, the values of potential at each impurity site might be considered as independent 
random variables. By homogeneity, they should have the same distribution. 

On the other hand, since the ground state wave function decays exponentially fast, the 
tunneling of a charge carrier between two impurity sites is controlled by the overlap of the 
two wave packets localized on each of the impurities. This overlap is of the order of e"'"/"^ 
if r is the distance between the two impurities. If the average distance between impurities 
is about lOa^, this term is very small (of the order of 10~^). This implies that the effective 
Hamiltonian describing the motion of a charge carrier through the impurity sites is made 
of two contributions: (i) the contribution of the random potential created by the ions in 
the crystal and (ii) the hopping term of the order of e"*"' "^ . Since this latter contribution 
is much smaller than the potential contribution, this effective Hamiltonian is given by 
an Anderson model in which the kinetic term is much smaller than the potential term. 
Hence, paradoxically, lightly doped semi-conductors, at very low temperature, correspond 
to a strong disordered Anderson model. In the present paper, the hopping term will be 
ignored, so that the one-particle Hamiltonian will be a pure on-site potential. 

2.4. Mott's Variable Range Hopping Transport. In 1968, Mott [32j proposed the 
following argument liable to explain the conductivity properties of Anderson insulators at 
low temperature (see also [32] )• Mott assumes that the solid is a d-dimensional Anderson 
insulator, namely the electrons are strongly localized. In particular, each energy level 
of the electron energy spectrum is associated with a position in the solid within a ball 
of diameter given by the localization length ^ (which can be taken of the order of a^). 
Moreover, he assumes that the density of states (DOS) np at the Fermi level is non 
vanishing but small enough so as to avoid overlap between electron states. In particular 
the mean distance between neighboring electron states is large compared to ^. The inverse 
of rip • ^'^ is a measure of the mean level spacing between states within a ball of diameter 
^. The temperature will be small compared with this scale, namely 

kBTnpi'' = — < 1, 

if kB denotes the Boltzmann constant. Then within a small error, all states with energy 
smaller than the Fermi level Ep are occupied, whereas the ones with higher energy are 
empty. The probability that a phonon of energy e is produced is proportional to e^^^B^ . 
This is true provided e ^ ksT. Such a phonon can be absorbed by an electron of the 
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Fermi sea localized at x G £, with energy e^ < Ep, and cause this electron to hop into 
a state localized at y E C with energy e^ > Ep. The probability for such an event to 
happen is controlled by the tunneling effect forcing the electron to move from x to y. 
Let r = |x — y| be the distance between such states, so that the tunneling probability be 
proportional to e"*"/^. Therefore, the probability P of transfer of an electron at distance 
r from its original location is proportional to 



P oc exp 



e r 



(2.3) 



By definition of the DOS, the product e ■ np ■ r'^ represents the average number of states 
in an energy interval of width e localized in a cube of size r. Mott then argued that the 
most likely value of the radius satisfies 



Up ■ e ■ r 



(2.4) 



The probability of jump P is then optimized over r and e. This can be done by giving r 
its minimum value compatible with eq. (2.4) and by maximizing w.r.t. e. Therefore the 



conductivity, which is proportional to the sum of contributions of all such jumps, will be 
controlled by the maximal value of P. Optimizing over e leads to 



P oc exp 



T 



i/{d+i) 



ksTo 



(rf+1) 



d+l 



d'^ np^° 



(2.5) 



Then, the phonon energy optimizing P is given by Eopt ~ 7"^/^+i ^ T, while the average 
distance of the jump is given by r/^ ^ (To/T)^/^+^ ^ 1. 

Clearly, the conductivity is proportional to the probability of transfer of electrons per 
unit time. Therefore, we expect the conductivity to be proportional to the same factor 
as a function of the temperature. 

For lightly doped 3D-semiconductors like silicon (see [39]), the localization length is given 
by the Bohr radius of the impurity, which is about lOOA. For a concentration of 10^^, 
the mean distance between impurities is lOOOA which is approximately 10^. Then, the 
typical width of the impurity band is about ImeV. Assuming the DOS to be flat on it 
gives Uf = c/W which leads to Tq ~ 1.1 x 10^ i^ ! This is huge indeed. At T = IK, 
this gives (Tq/T)^/^ ~ 18, namely (i) the electron hops at about 18^ ~ 2rat,, (ii) the 
conductivity is multiplied by a factor e"^^ ~ 1.2 x 10~^ due to this mechanism ! It turns 
out that this is exactly what happens in the quantum Hall effect (QHE) |3H [5]: the Mott 
variable range hopping controls the fluctuation of the plateaus, leading to the amazing 
accuracy of this experiment. 

Mott's law has been well documented in the experimental physics literature. A review 
of the results obtained in the seventies on various semi-conductors can be found in [22] 
(see also [53]) and some results have been reported in [33]. In addition a large part of the 
book by Efros and Shklovskii [39j is dedicated to this effect (see Chapter 9, in particular 
and references therein) . 
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The previous version of Mott's argument is based on a critical assumption summarized 



by eq. (2.4). This part of the argument has been the focus of attention of several works 
in the early seventies [Ml [ISl El I2SI EZl ESI E] and it is the main topic of a large fraction 
of the book of Shklovskii and Efros [39]. Following the description provided by Miller 
and Abrahams [31], the electron conduction is seen as a random resistor network. In 
the early sixties, Ziman [5l] suggested that percolation theory should be a key technique 



to investigate hopping transport. Since then it became indeed increasingly clear that 



percolation gives the right argument to justify eq. (2.4) and Mott's prediction 



2.5. Description of the Model. Based upon the argument of Mott, it becomes possible 
to propose a model. This section will be dedicated to its heuristic description, leaving the 
rigorous definitions for the Sections that follow. All along the present paper, only n-type 
doped semiconductors will be considered. A particle-hole symmetry permits to consider 
the p-type case. 

Let L be the lattice TL^. For any x G £ let s^; G {0, 1} be a random variable with s^, = 1 
if and only if an electron state is available at x. Given a family of random variables 
■s = (sa;)xe/: each taking values in {0, 1}, let £(s) be the random subset of L containing 
all sites x where s^ = 1- Thus C{s) is the set of sites in the semiconductor on which an 
impurity electron state is available. In the tradition of Solid State Physics, and only for the 
heuristic description of the model, the total number of sites available will be considered 
as finite. The electron will be described in the second quantization picture, (see [16]). 
through a pair of Fermion creation-annihilation operators [a]., ax), where x G C{s) is 
a lattice site, thus obeying the canonical anti-commutation rules which we abreviat by 
CAR, 

a^ciy + ttyttx = , a'^ = a]. , a^dy + ayttx = S^yl ■ (CAR) (2.6) 

Here A* denotes the adjoint of the operator A. In this work, for simplicity, only spineless 
electrons are considered, because the model proposed will not couple the spin to the mo- 
tion. The quantum dynamics will be made of two parts: (i) the coherent part, describing 
the electron motion in absence of dissipation, and (ii) the dissipative part, taking into 
account the interaction with other degrees of freedom, especially the phonons. 

Since the dynamic concerns the electron gas, the system will be considered in the local 
equilibrium approximation ^ . Namely the electron gas is seen as the union of mesoscopic 
cells with size large enough so that within a mesoscopic time each such cell has the time 
to return to equilibrium. On the other hand the size of the cells is small compared to 
the size of the sample under consideration. In each mesoscopic cell, the electron gas is 
in an equilibrium state, namely defined by the usual thermodynamic parameters, the 
temperature T (or rather the inverse temperature (3 = l/fc^T) and the chemical potential 
/i. Such mesoscopic cells are opened, namely they allow both energy and electrons to be 
exchanged from cell to cell. Therefore, only the average of the energy and of the particle 
number is fixed. In addition, if the system is put out of equilibrium, then both T and /i 
may vary slowly in space and time, so that the relative variation in each mesoscopic cell 
is negligible. 
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The coherent part is given by a Hamiltonian and the observables evolve according to the 
Heisenberg equation. But, because of the local equilibrium approximation and because 
each mesoscopic cell is opened, the Hamiltonian will represent the free energy of the gas 
in the mesoscopic cell. Therefore it has the form F = H — fiN where H represents 
the mechanical energy of the gas, while A^ is the number of electrons (in the mesoscopic 
cell). The mechanical part H of this Hamiltonian will be reduced to the potential energy 
of the electrons on the impurity sites, ignoring the hopping term between impurities, 
since the tunneling effect is so small for lightly doped semiconductors. The potential will 
be represented by a family e = (ex)xG£(s) of independent identically distributed random 
variables, where e^. belongs to the impurity band A which is a compact interval. For the 
sake of the present model, the distribution will be assumed to be uniform in A and it 
shall be independent of the random variables s = {sx)xgc- 

The family of all random variables s = (sx) and e = (e^^) will be simply denoted by a 
random variable u = (s, e). The corresponding probability space consisting of the compact 
set of all such families will be denoted by Q and the corresponding probability measure 
on the Borel sets of Q by P. Many operators such as the free energy depend on the 
randomness u, e.g. the free energy becomes 

F{uj)= J2 (ex-/i)4«x (2.7) 

x'6£(s) 

Since the electron-phonon interaction is weak, the dissipative part will be given in the 
Markov approximation, in which all degrees of freedom other than the electrons are in- 
tegrated out, while the time scale for this interaction to rearrange the electron state is 
considered as negligible. Namely, the dissipative part of the dynamics will be described 
by the generator of a Markov semigroup, which, thanks to the Theorem by Lindblad f28] 
is given by a Lindbladian of the form 



2^(^) = E (l(L^L,A + AL*L,) - L*AL^ 



(2.8) 



The description of the model will precisely consist in proposing an expression for the Lj's. 
In order to implement the Mott scheme, these operators should describe the jump of an 
electron form an impurity site x G C{s) to an impurity site y G C{s). The corresponding 
jump operator, denoted by L^^y, should be proportional to aj^Oa,, since the later annihi- 
lates an electron at x and creates one at y. If one of the sites, x or y, is not in the random 
lattice jC(s), then there can't be a jump from x to y. Hence we set L^^y equal to zero in 
this case. Hence, 

Lx^y{uj) = SxSy a/TV^ alax . (2.9) 

where T^-^y is the probability rate at which the jump arises. Following the Mott argu- 
ment, this jump probability rate needs to take into account the probability of a phonon 
being absorbed or created by the electrons. The absorption process is dominated by the 
probability of a phonon of energy e to be created by the thermal bath. It is given by 
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the Boltzmann factor e"'^'^ with a good approximatioKQ Then the electron located at x 
absorbing this energy must find an available site y with energy e^ = e + e^ > e^- On 
the other hand, if an electron at x is in an excited state, it might decrease its energy 
by spontaneously emitting a phonon or energy e > 0, provided it finds a site y at which 
the energy available is e^ = e^; — e < e^^. The spontaneous emission does not require any 
Boltzmann factor. Hence the probability for absorption and emission is proportional to 
Q-P{<^y-<^x) jf g^ > g^ g^]^^ ^Q 1 if ej; < ex- Thus a unified formula for the probability of 

absorption or emission is e"^*^"^^"*^""-* where e"*" denotes the positive part of the real number 
e. This difference between emission and absorption leads to the quotient p^^^ = e~^^'^y~'^'^^ 
which is also known as detailed balance condition. 

In addition, in both cases, the electron must jump form x to y through a tunneling 
effect, namely decreasing the probability by a factor proportional to e~''^'~^'/'', if r is the 
localization length. Note that in lightly doped semiconductors, r is of the order of the 
Bohr radius a^, since the average distance between impurities corresponds approximately 
to lOttB- This proportionality factor will be normalized by dividing it by Z, where 

e-l"^l/^ (2.10) 

This leads to the expression of the jump rate proposed by Mott 

— \x — 111 1^ 

r,^,H = To '—^— e-P^'y-'^^^ , x,i/ G C{s) (2.11) 

which gives the detailed balance condition 



E 



r 



e-P^^y-^^) . (2.12) 



Remark 2.1. In the 'present paper we assume that Tq is a constant. The present results 
however remain valid if To is a function of ex, ey and /3 satisfying the following properties: 



To{ex, ey, j3) = To{ey, ex, /?), (hence the detailed balance condition (2.12) will remain valid), 
and inf{ro(ea:., ey, /3) : e^, e^^ G A, /3 G (0, oo)} > 0. 

The factor Fq is a parameter with the dimension of the inverse of a time, fixing the order 
of magnitude of the effect. We let 2)^*" denote the Lindblad operator obtained from this 
family of jump operators, 

S''"(^) = J2 il {LUyLx^y"^ + ALUyL.^y) ' LUy^Lx^y) ■ (2-13) 

However, this part of the model is insufficient to describe the return to equilibrium, 
because S)*'*"' leaves the number operator N invariant, as can be checked easily. Therefore 
it does not take into account the possibility for an electron to jump out or to jump in the 
system (thermal bath). An electron can jump out of the system in two ways: (i) either 



It should actually be given by the Bose-Einstein distribution [e^"^ — 1) ^ . But if /3e ^ 1 it follows 
that (e'^" - l)-i « e-<^^ 



12 GEORGE ANDROULAKIS, JEAN BELLISSARD, AND CHRISTIAN SADEL 

its energy becomes too large or too small to stay within the impurity band A, or (ii) the 
electron is kicked out of the mesoscopic cell under consideration. Similarly the opposite 
processes arises to allow an electron to jump in the system. A standard way to take this 
process into account is to create an extra site ^ called the cemetery and to describe these 
processes as a simple jump ^ — )■ a; or x — )■ ^. In the first case, an electron is created at x, 
while in the other it is annihilated at x. The cemetery really describes the thermal bath 
and it is only natural to interpret the chemical potential fi as the energy associated with 
this new site. Hence, the corresponding jump operators will be given by 

L^^^(w) = s^ (r^_,^)^''^ a^ , L^_,^(w) = s^ (r^^^)^/^ 4. . (2.14) 

A convenient way to define the jump rates T^^^ is to imitate what was done earlier and to 
consider the sites * similar as the sites occupied by electrons, associated with the energy 
fi, leading to 

T^^^ico) = r. e-''('^--)^ , T^^^iu) = r, e-^^--'^)^ , ^ = e^^-"'^) (2.15) 

where x G C{s). Here again, F^ > is a jump rate probability fixing the time scale for 
the cemetery process. 

Remark 2.2. In the present paper T^ is a constant. The present results however remain 
valid ifT^, is a function of e^, and (3 satisfying inf{r^(ea;,/3) : e^. G A,/3 G (0, 00)} > 0. 
See also the related Remarks \2.1\ and \5.9[ 

The new jump operators describing the exchange between the thermal bath and the 
electron gas gives rise to a Lindblad operator denoted by 2)*. However, due to the anti- 
commutation rules, the Lindbad operator, acting on an observable A, takes the form 



KiA) 



Y, (^ (LU^L^^^A + ALU^L,,^,) - i-l)"- LU^AL,^.^ + 



x<=C{s) 



J2 (l (LU^L^^^A + ALU^L,^.,) - {-lY--LU,AL,^}i (2.16) 

where c?^ is the degree of A given by the natural Z2 grading of the CAR algebra 2l(s). 
More details will be given below. Hence, the total dynamics is described by an operator 
£ acting on the set of observables by 

£.(A) = i[F{u:),A] + D^-(A) + ^1{A) . (2.17) 

What is left for the mathematician, is to make sure that this description does not produce 
any hidden effect that could be related with the infinite volume limit. For indeed, the 
mesoscopic cells have an undefined size, only the order of magnitude of the size is fixed. 
In addition, the random character of the set C{s) of the impurity sites and of the e^j's 
must be included in the description to make sure that the results obtained are almost 
surely independent of the configuration of the impurities. These two requirements are the 
very reason why the mathematical formalism is so demanding. 
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3. The Coherent Evolution 

In this section the coherent evolution of the quantum motion is studied and the equihbrium 
state of the unperturbed system is obtained. The coherent evolution is a one parameter 
group automorphism on the algebra of the observables. 

3.1. Observables. In this subsection the CAR algebra of the observables is studied. 
Recall that a CAR algebra is a C*-algebra generated by a countable number of creation 



and annihilation operators which satisfy the rules stated in eq. (2.6). It is well-known, and 
elementary to show, that the complex algebra generated by the pair ax-ia)^ is isomorphic 
to the set M2(C) of 2 x 2 matrices by using the analogy 





1 



a;^^ 



1 




If a finite family A C C{s) of sites is considered instead, the Z2-graded algebra generated 
by the family {a^, aj. ; a; G A} is isomorphic to the tensor product M2(C)®^ but the previ- 
ous analogy ought to be modified in order to insure that ax and ay anti-commute. One way 
to describe such an isomorphism is given by the so-called Jordan- Wigner transformation: 
since A is finite let its elements be numbered from 1, 2, ■ ■ ■ ,m = |A|. Then 



a-k ^ o^ 



0-3 (8)a (g) I2 (8) 



fe-i 



^k^ZL 



(Ts (X>cr 



I2® 



(3.1) 



fe-1 



where 



c^a 



1 
-1 



As can be seen, the Jordan- Wigner transformation is not canonical in that it requires the 
choice of an order on A. However it can sometimes be convenient for practical uses. 

A canonical description of the observable algebra has been given in the past and the main 
references are [9l[T0]. Let this construction be summarized here. As shown before, there is 
a non-canonical isomorphism between 21a(s) and M2{C)'^-^'^^^'^\ so that ^a{s) has dimen- 
sion 2l^'^^'^'')L If A C A' there is a canonical isometric embedding ^aa' : ^a{s) ^ ^a'{s) 



(in the Jordan- Wigner representation (3.1) this embedding consists simply in adding to 
the Ofc's one more factor ®l2 at the end of the chain). The algebras 21a(s) are called local, 
while the inductive limit 

2l(s)=lim(2lA(s),2A,A') , (3.2) 

is equal to the quasi-local observable algebra. It should be remarked at this point that 
2l(s) is random as it depends upon the random variable s = {sx)xec- However, since 
s G {0,1}^ = H and since S can be considered as a compact space, this family of algebras 
can be seen as a continuous field of C*-algebras [TTl STJ HHl HH]. 

Recall that a C*-algebra algebra 21 is Z2 graded when there exists a ^-automorphism 
a on 21 which satisfies a^ = 1^, (in other words, when the group Z2 acts on 21). The 
^-automorphism a is called grading. In the case of the CAR algebra, the grading is the 
canonical ^-automorphism defined by a{ax) = —ax for all x. If 21 is a Z2 graded C*-algebra 
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with grading a then the elements A G 21 with cr(A) = A are called even, and the elements 
A e 21 with (j{A) = —A are called odd. The degree of every even element A is defined to 
be equal to zero, (denoted by (Ia = 0), and the degree of every odd element A is defined 
to be equal to one, (denoted by cIa = 1). In the above Jordan- Wigner representation, 
the operators of degree zero are represented by diagonal matrices, while the operators of 
degree one are given by the off diagonal ones. 

Proposition 3.1. The family {2l(,s) ; s G S = {0, 1}^} can be endowed with the structure 
of a continuous field of'L2-graded C* -algebras. 

Proof: Recall first ^^7\ that a section of the field (2t(s))^g- is an element of the Cartesian 
product rises 2t(s). To define a structure of a continuous field it is necessary to define first 
a set "J of reference sections with the property that (i) for each s G S the set {^(s) ; ^ G 5"} 
is dense in 2t(s) and (ii) if ,^ G 3", the map s G S i-t- ||^(s)|| G M+ is continuous. Then a 
continuous section is an element rj G Hse^ ^('^) such that for any s G H and any e > there 
is a neighborhood U of s in S and a reference section ^ G 5" such that \\^{s') — ri{s')\\ < e 
for s' G U. 

In the present situation it ought to be remarked that the product topology on S is defined 
through the set of open sets 'Ua(s) = {s' G S ; s'^ = s^ ,Wx E A}. Hence s' is close to s if 
there is some finite subset A C £ such that s'^ = s^ ior x E A. Then let a^ be the section 
defined by a^i^s) = s^a^ G 2t(s). Let 5" be the *-algebra defined by the a^^s. Namely 
^ G 5" if and only if there is A C £ finite, an integer M G N and continuous functions 
{Xtu-'-bm ^ ^(^)} such that 

M 

^(^) = E E E ^" i',i (^) ^-1 ■ ■ ■ ^-^ < ■ ■ ■ 41 (3.3) 

m=0 Xl,...,Xm&^ ttl vflm 

where the jjj's G {., f} label the a^^, a^^s. Clearly, if A C A' and if s' G 1Xa'(s) the previous 
expression for ^(s) does not change as s' replace s, so that s G S i— > ||^(s)|| G IR+ is 
continuous. Moreover, the set {^(s) ; ^ G 5"} generate 2t(s) as a C*-algebra. Hence 3" 
satisfy the conditions. □ 

Proposition 3.2. The continuous field ((2l(s)se=, 5") is covariant by the translation group, 
namely, for any a G Z'^, there is a 7j2-graded *-isomorphism rja : 2l(s) i— )■ 21(t"s) leaving 
3" invariant. 

Proof: It is elementary to check that the ^-isomorphism rja is defined by ria{ax){s) = 
Sxttx+a exists and does the job since {T"'s)x+a = Sx- Moreover it is easy to check that 
Va° Vb = Va+b, that Tjo = id and that rja commutes with the grading. □ 

The continuous field ((2t(s))ses, 3^) contains an important sub-field made of commutative 
C*-algebras. For given s G S, let £(s) be the closed sub-algebra of 2t(s) generated by the 
family {sxUx ; x G £} where rix = a\.ax. It is easy to see that (t(s) is commutative, it 
contains only elements of degree zero, it is generated by projections and generates also a 
continuous Z'^-covariant field of C*-algebras. 
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Two new ingredients will be needed in the rest of the paper. First it will be convenient 
to consider the field 21 as a field over the compact space Q instead. The only difference 



will be that, in the eq. (3.3), the coefficients will be allowed to depend continuously 
upon u = (s,e). When necessary this field will be denoted by 2l(a;), and similarly the 
corresponding commutative sub-field will be denoted by €{u). The following can be found 
in [Hj. 

Corollary 3.3. The field 21 can be continued as a continuous field 21 = (2l(a;))^gj^ of C*- 
algebras which is covariant with respect to the translation group. Similarly € = {€{u))^^^ 
will denote the abelian sub- field extending {€{s))^^.= . 

The other ingredient is the use of groupoids to describe the Z'^ action [HI EH]- Here 
the group Z'^ acts on the space Q by t°-{uj) = T°-{s,e) = (t''(s), T''(e)) where a G Z"*, 
{t'^{s))x = Sx-a and (T"(e))a; = ex-a- Consequently, the crossed product Fq = fi x Z'^ 
is a groupoid described as follows: (i) elements are pairs 7 = (w, a) G i7 x Z'^, (ii) Q 
is the set of objects or of units, (iii) each element has a range and a source in Q, here 
r{u,a) = u , s{u,a) = t'^-oj, (iv) two elements 7 and 7' are composable if 5(7) = ^"(7'), 
namely if 7 = {u, a) then 7' = {t~°'u, a') and there is a composition law 707' which 
is here given by (1^,0 + a'), (v) the elements of the form (c<;,0) are units and can be 
identified with the points in Q, (vi) each element 7 admits an inverse with exchange of 
range and source, namely here 7"^ = (T~"a;, —a). The groupoid Fq will be endowed with 
the product topology and it is elementary to check that all groupoid maps defined above 
are continuous. 

To express the covariance of the field 21 = (2t(a;))^gf^, it is convenient to see the translation 
Tja as a function of the groupoid variables instead, namely if 7 = {u,a), then 7 can be 
seen as sending its source 5(7) = T~'^a; into its range r{'-f) = u, so that 

V{u,a) : 2l(T-'^a;) ^ 2t(u;) . 

Definition 3.4. A field (^'lj)^^^^ of *-homomorphisms from 2t(aj) into itself will be called 
covariant if the following diagram is commutative 

2t(T-'^u;) "-^^ 2l(w) 

2t(T-'^u;) "-^^ 2l(w) 

This field will be called continuous whenever it transforms every continuous section of 
Qi{u) into a continuous section. 

3.2. The coherent dynamics. Recall that the coherent part of the quantum motion is 
a group automorphism acting on the CAR algebra 21. Let A C £ be a finite set. Then 
the particle number A^a, the Hamiltonian if a and the free energy Fa within A are defined 
as follows: 

^a(^) = ^Sxnx, Hf^{u) = ^ s^ex n^ , (3.4) 

x£A x£A 
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Fa(w) = Ha{uj) - fiNA^u) = ^ Sx{e^ - fi)nx , (3.5) 

xeA 

where 

Hx = a\.ax and w = (s, e) G fi . 

The finite volume coherent dynamics is generated by the free energy operator 

aJ"'^^(A) = e**(^A-/.AfA) ^ ^-^t(H^-^.N^) ^ ^ ^ ^1^(0;) , t g M . (3.6) 

Usual arguments [TU] , which become trivial in the present situation, show that the infinite 
volume limit exists, namely 

af\A)=\ima'f'^\A), A G 2tAo(s) , t G M. (3.7) 

Indeed, it is elementary to show that 

at\sxa^x) = e^'^^^'^-^^sxal , at'\sxax) = e-^'^^^'^-^hxax , (3.8) 

since it is already true for af'' as soon as A 3 x. Hence, of can be computed on any 
monomial, thus on any polynomial, in the creation-annihilation operator. Therefore af 
is defined anywhere in 2t(s). It also follows from these formulas that s^fix is invariant 
by the dynamics, a fact which comes from ignoring the hopping terms due to possible 
tunneling between impurity sites. As a consequence, the elements of €. are left invariant 
by the dynamics. To summarize, thanks to the Definition 3.4, the following holds (the 
proof is left to the reader) . 

Proposition 3.5. The field of equilibrium dynamics a'-'^^ = («( ) induced by the 
Hamiltonian (3.4) on each 2t(a;) is well defined, iZ-covariant and continuous. Moreover 



its restriction to the sub-field €. is trivial. 

The generator of the coherent evolution af is given by the commutator with the free 



energy operator defined in eq. (3.5), i.e. 

dat\A) 



5M) 



dt 



t=o ^t/: 



\iim[FA{uj),A] for Ae%oc- (3.9) 



3.3. Equilibrium State. Once the dynamics is defined, the next step is to find the 
possible equilibrium states. There are two ways to do that. The first is to consider the 
finite volume approximations and establish the quantum version of the DLR equations. 
The other one, valid in the infinite volume limit, consists in implementing the Kubo- 
Martin-Schwinger conditions (KMS) which were introduced in [26l[30] and further studied 
in [21I HH [ini HE] . In the present situation, due to the extreme simplicity of the dynamics, 
both approaches can be used and lead to the same explicit result. The KMS condition 
will be used here. 
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Definition 3.6. Let p he a state on aC*-algebraQi, and let a = {at)t£R be a one-parameter 
group of *- automorphisms o/2l. Then p will be called P-KMS with respect to a if 
(i) it is invariant by the dynamics, namely po at = p for all t e M, 
(a) if A,B e 21, then 

p{AB)= p{a.^p{B)A) (KMS condition) (3.10) 

Note that (i) is a trivial consequence of (ii) if 21 is unital, but it is convenient to define 
this condition in this way. In the present situation the following holds. 

Proposition 3.7. For ^-almost every u = (s,e) G Q, there is a unique /3-KMS state p^ 
on 2l(co') with respect to the dynamic a*^'^\ It is defined by 

A,. G2t|,|(s) (3.11) 

(3.12) 

(3.13) 

Moreover, the field of such states is continuous and covariant with respect to translations, 
namely p^ o r]^^^ = p^-a^. 



11-4. = 


Puj{A^) , 


KxeA / 


xeA 


Pa;(Sx4) = 


= 0, 


puisxrix) = 


1 


"^■l + e^(^--^) ■ 



Proof: (i) In order to prove eq. (3.12) and (3.13), it is sufficient to assume that s^ 



Then, each element of the elementary algebra ^{x} can be written as a linear combination 



of l,ax,al. and n^ = al.ax- Thanks to eq. (3.8) and since a KMS-state is time invariant, 
it follows that PujicLx) = Puiio!t{ax)) = e~**(''^~^^p^(aa;) for all t G M. Since the distribution 
of Sx is absolutely continuous, ex ^ p P-almost surely. Hence p^{ax) = 0. In much the 
same way PujiO'l) = 0. Similarly, ii A,B ^ ^a^uj) where x ^ A, the same argument shows 
that pi^{AaxB) = pi^{Aal.B) = 0. Consequently if A G 2l(a;) \ ^{uj) it follows by induction 
on A and by density, that Puj{A) = 0. Hence it is sufficient to reduce the analysis on the 
commutative sub-algebra £(a;). 

(ii) Since Ux commutes to any element of €a{u), it is sufficient to consider an expression 
of the form Puj{Anx) = puj{Aal.ax) where A is a local observable in €.{ijj) with support not 
meeting x. Using the KMS-condition, this gives 

Puj{Anx) = Puj{a-^fj{ax)Aal) = e"'^(^""^Va;(«x^4) • 

Since the support of A does not meet x and A G €{uj), it follows that A commutes with 
Gx so that puj^axAal) = p^^Aaxol) = p^{A{l — Ux)). Hence this gives 

o (An)- ^"^"^^ 

Eliminating points of the support of A one after another leads to the formulas above. 

(iii) It is obvious that the field p = (ptj)^gn ^^ states is continuous. For indeed it is enough 
to show that the map u E Q ^-^ P(^(^(a;)) G C is continuous for any continuous section 
of the field 21. By definition of continuous sections, it is sufficient to chose (^ G 3". But 
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this is exactly choosing ^(w) as a polynomial in the creation-annihilation operators. The 
formulas (3.11), (3.12) and (3.13) show immediately the continuity with respect to u. 



(iv) Similarly, the covariance with respect to translation follows immediately from the 
formulas (IS^, KV2}^ and (Islsl). □ 



Proposition 3.8. For all u E Q the state p^i is faithful. 



Proof: First the state p^ is defined by the formulas (3.11), (3.12) and (3.13). It is easy 



to check that such result can also be obtained from (omitting the reference to u) 



PiA) 



Tr (e-'^^^A) 



Ae^A 



where Z{A) is a normalization constant and Fa = Ha — ^Na. In particular, the restriction 
of p to SIa is faithful. Moreover, there is a conditional expectation Ea : 21 t— )■ SIa defined 
by 

TrAvle-^^'^'vM) 



Ea(A) = Ea',a(A) 



if A G 21. 



It is elementary to check that A > ^ Ea(A) > 0, that Ea(1) = 1, so that ||Ea(A)|| < 
1 1 A 1 1 for A G 21. Moreover, it is easy to check that p o Ea = p for all finite A C £. 
Let now A G 21 be positive and such that p{A) = 0. Then p{A) = p(Ea(A)) = 0. 
Since p is faithful on 2tA it follows that Ea(A) = for all finite A's. Let e > 0, there 
is Ae C £ finite and A^ G 2tAg such that \\A — A^W < e/2. Hence, whenever A D A^, 



|Ea(v4 - Ae 



\As\\ < e/2. Therefore ||y4|| < e for any e implying A = 0. 



D 



3.4. The GNS Representation. The Gelfand-Naimark-Segal construction (GNS) is a 
fundamental tool in the study of C*-algebras. Let 21 be a C*-algebra, for convenience it 
will be assumed to be unital. Let p be a state on 21. Then a Hilbertian inner product can 
be defined through {A\B) = p{A*B) for A,B eOi. The set 9T of elements A G 21 for which 
p{A*A) = is a closed left 2t-module. The quotient space 2t/9T inherits the structure of 
a separated pre-Hilbert space. By completion it gives a Hilbert space Ti = L'^{^,p) and 
a canonical linear map ( : A E'Ql^-^ C(^) ^ ^ such that 



{aAmB)) = p{A*B) 



(3.14) 



If 21 is abelian, then, by Gelfand's theorem, it is isomorphic to the set of continuous 
functions on some compact space X, unique up to homeomorphism, called the spectrum 
of 21. Then a state is just a probability measure on X and T-L = L'^{X, p). 

% inherits the structure of a left 2l-module so that the map 7r(y4) : C{.B) eH ^-^ (i^B) G 
T-t is well defined and extends to a representation of 21 in 7{. The vector ^ = ({!) is cyclic 
since C(^) = ^(^)'C so that H is the closure of the set of 7r{A)^ as A varies in 21. The 
weak closure of 7r(2t) is a von Neumann algebra denoted by L°°(2t, p) = Tl. Clearly, p 
extends as a normal state on VJl since 



p{A) = {^\n{A)0 . 



(3.15) 
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If 21 is abelian with spectrum X, the representation is given by point wise multiphcation 
in L^ {X, p) . The vector ^ is the constant function equal to one and the von Neumann 
algebra is the space L°°{X,p) of essentially bounded p- measurable functions on X. 

If a is a ^-automorphism of 21 leaving p invariant, it defines a unitary operator Ua '■ 
C{A) ^ ({a{A)) G n such that 

UaniA)U-^ = n{aiA)) , A G 21 , U^^ = ^ . (3.16) 

The Tomita-Takesaki theory [H] is based upon the densely defined conjugate linear opera- 
tor S defined by S({A) = C(^*)- The main result of Tomita and Takesaki is that 5" is well 
defined, closable and that, if S denotes also the closure, A = S*S is a positive self-adjoint 
operator on 7i called the modular operator. The polar decomposition S = JA^/^, defines 
a conjugate hnear involution J such that JWIJ coincides with the commutant of 971 on 
T-L. It has been shown [H] that there is a group of *-automorphism 6 on DJl satisfying a 
KMS-condition, namely 

OtiM) = A**MA-** , p(MiM2) = p(e„,(M2)Mi) , M, Mi, M2 G 971 . 

In the commutative case the modular operator is trivial while JA coincides with the 
complex conjugate of A. In the early seventies, Araki [1] and Connes [13] realized that 
the set of C(^) with A G 21, A > generates a cone 71+ defined by 



H+ = {AV%(A)e; AG2t, A>0}, Jn+ = U+. 

Connes characterized such a positive cone in L^(2t, p) as being self-dual, homogeneous 
and oriented. It is the non-commutative analog of the set of positive square integrable 
functions. 

In the present situation the general theory gives the following result. 

Proposition 3.9. Let 21 = (21(0;)) „ he the continuous field of C*- algebras defined in 
Section 3.1 Let a = (a'-'^^) be the field of dynamics defined in Section 3.2. Then the 
GNS-construction leads to a continuous field l-i = {'Hu})u;en of Hilbert spaces, a continuous 
section ^ = (^a;)wen of unit vectors, with a U^ action, namely a unitary representation of 
the groupoid Fq, given by unitary maps K;,a : 'K-^-'^w '— ^ '^^ satisfying 

The 'Ij2-grading in 21 is represented by a continuous covariant field G = (G'(^)a;en of oper- 
ators satisfying G^^ = G*^ = G~^ and the covariance condition 

The field 21 is represented by a continuous field n = {nca)ujen of representations, for which 
^ is a field of cyclic vectors, and satisfying the covariance condition 

VL,aVrT-^(A)K.J = 7r^(r/^,a(A)) , A G 2t(T-'^a;) . 

The dynamic gives rise to a strongly continuous group of unitary operators on 71^ with 
generator F^. The latter defines a continuous covariant field of self-adjoint operators 
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satisfying 

The family of Tomita- Takesaki modular operators define the covariant continuous field of 
KMS-dynamics as follows, 

where jj denotes {-,1} in the exponent of the a's, while it corresponds to { — 1,+1} respec- 
tively in the exponential. The Araki-Connes cones 7/1^,+ gives also a continuous covariant 
field of self-dual, homogeneous, oriented cones generated by vectors of the form 

Since the proof is straightforward it will be left to the reader. 

Remark 3.10. The notation F^ for the generator of the dynamics in Tiuj is justified, 
because it corresponds to the free energy and can be seen as an infinite volume limit of 

4. DissiPATivE Dynamics 



The general dissipative operator 1)^ = 2)^*" + 2)* defined in Subsection 2.5 is the main 



object of study of this section. 2) can be considered as a field of operators {1)ui)ui€n on 



the field of C*-algebras (2t(a;))^gn. Theorem 4.9 is the main result of the section where 
the Friedrich extension Theorem is used to prove that the closure of the operator 2)(^ 
is a positive self-adjoint operator on the Hilbert space of the GNS representation of the 
algebra of the observables using the equilibrium state. 

4.1. Complete Positivity. If 21 is a unital C*-algebra then a map r^ : 21 — )■ 21 is positive 
whenever A > ^ vi^) ^ 0- Then t] extends to 21 ® M„(C) by ?7„ = r^ ® id. Then r] 
is called completely positive, (a term which was introduced in |13]), if ?7„ is positive for 
all n's. rj is called normalized if ?7(l2t) = I2I) whenever I21 denotes the unit of 21. By 
CP(2l) we will denote the set of completely positive maps and by CPi(2t) the subset of 
normalized CP-maps. Examples of CP-maps are: 

(1) Any *-homomorphism is CP. It is normalized if it sends 1 to itself. 

(2) If L G 21 then the map A H- L*AL is CP. It is normalized if and only if L is an 
isometry. 

(3) If $ is CP and $(1) is invertible, then the map <I>'(A) = $(l)-i/2$(^)$(i)-i/2 jg 
CP and normalized. 

(4) Any convex combination of CP-maps is CP and the same holds for CPi. 

(5) The composition of two CP-maps is CP and the same holds for CPi. 

(6) If ($n)neN is a sequence of CP-maps such that $(A) = lim„_j.oo *&n.(^) exists for 
all A G 21, then $ is CP, the same holds for CPi. 

(7) If L G 21 let ^l(v4) = L*A + AL. Then e-**^(A) = e-*^*Ae-*^ is CP for t > 0. 

As a result we get the following. 
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Proposition 4.1. Let Qi be a 1j2-graded C* -algebra with grading automorphism a. Let 
\1> : 21 ^-)■ 21 6e a linear map of the form 

N 

^(^) = E { 2 (^^^'^ + "^^^'^ ~ i^ir^^'-LlAL,] , (4.1) 

where Li & 01 for all i 's and all of them are either even or odd elements. Then the map 
$i = e~** commutes with the grading and is CPi for t > 0. 

Proof: Since \1/ is a bounded operator, the exponential does exist. Moreover ^1^(1) = 
since the degree of 1 is zero, so that $i(l) = 1 for all t's. Since \I/ is a finite sum of 
operators, the Trotter product formula [501 El] will prove that $t G CP as soon as each 
pieces of the sum gives a CP-map. Since the grading is given by a ^-automorphism, 
the map A i— )■ {—lY^i'^'^L*ALi can be seen as the composition of two CP-maps namely 
a I—)- a'^{A) (with d the degree of Lj) and A i— )■ L*ALi. Exponentiation is given by a limit 
of polynomials in these maps with positive coefficients, namely it is CP. The other maps 
have the form A i— )■ L*LiA + AL*Li, the exponential of which is CP as well (see the last 
example in the list above). □ 



The proposition above justifies the form of the Lindbladian in eq. (2.8 ) for the generator of 
a Markov semigroup acting on a finite C*-algebra. It also gives the extension to Z2-graded 
C*-algebras . 

Let now p be a grading invariant state on 21. Then the GNS-construction gives a Hilbert 
space "H = L^(2t, p), a representation of 21 and a cyclic vector ^o- Moreover, the state 
p is cr-invariant, so that a defines a unitary operator G on "H, called the degree, such 
that G^o = ^0- Then clearly G7r(a)G~^ = n{a{a)). In addition, since a is an involution, 
G^ = 1^, G* = G. Let A, J be the modular operator and the modular involution (see 



Section 3.4) and let "H^ be the Araki-Connes homogeneous self-dual oriented cone in "H. 

If 21 = M„(C) and if p is the normalized trace tr„ then the Hilbert space "H is the 
space L'^{Mn) of Hilbert-Schmidt operators on C". It can be seen as the set of families 
{^ij)i<i,j<n with inner product {x\y) = ii/n)'^--Xijyij. The corresponding modular 
operator is trivial while the conjugacy is defined by jn{x)ij = xji. The positive cone is the 
set HS+(?7.) of positive nxn Hilbert-Schmidt matrices. Consequently, M„(2t) = 2t®M„(C) 
can be endowed with the states p ® tr„, leading to the Hilbert space l-L ® L^{Mn), made 
of families (xjj)i<jj<„ with Xij G Ti. The inner product is {x\y)n = {^/n)J2ij{^ij\yij)- 
The modular operator is then A^ = A ® 1„, with modular conjugacy J (g) j„, namely 
Jn{x)ij = J(xji). The corresponding positive cone will be denoted by 'Hn+- 

Definition 4.2. A bounded linear map F on Ti is positivity preserving if FiTi^) C 71+. 
It will be called completely positive if Fn = F ® 1„ is positivity preserving on Tin for any 
n. 

It follows from this definition that if (-Fm)meN is a sequence of completely positive maps 
on T-L converging weakly to F, then F is completely positive as well. 
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4.2. Markov Semi-Groups. One important property of tlie Lindbladian is the analog 
of the Leibniz formula for the second derivative, namely (fg)" — f'g — fg" = 'if g' ■ The 



following formula shows that a Lindbladian as defined in eq. (4.1) behaves like the analog 
of —A if A is a Laplacian. For the statement of this proposition, if 21 is a Z2 = Z/2Z- 
graded CAR algebra, define the graded commutator by 

[A, B\g = AB- {-if^'^^BA ioT A,B eOl. (4.2) 

A graded derivation of degree d G Z2, is a linear operator on the CAR-algebra such that 

6{AB)=6{A)B + {-iy-'^^A6{B). (4.3) 

In particular, if X is an element of the CAR-algebra, the map 6x '■ A ^^ [X, A]g is a 
graded derivation with degree dx- 

Proposition 4.3. Let £ be the Lindblad operator given by £^{A) = t[F,A] + \& where 



F = F* and \E' is given as in eq. (4-1)- Then for any pair A, B E 01 

(t) z{A*) = {ii{A)r 



N 



(11) Z{A*B) - A*Z{B) - £(A*)5 = - T.i=i[Li,A]*g[U,B]g (Leibniz formula). 

Proof: (i) Since F = F*, the first claim is obvious by inspection. 

(ii) The map 6 : A E ^ ^-^ t[F,A] is a ^-derivation and therefore it satisfies the Leibniz 
formula, S{A*B) = 6{A)*B + A*6{B). Hence this part does not contribute to the r.h.s.. It 
is enough then to consider the case A^ = 1. Let J denote the left hand side of the Leibniz 
formula for A^ = 1. Then 

2J = L*LA*B + A*BL*L - A*L*LB - A*BL*L - L*LA*B - A*L*LB 

_2(_i)dL(dA+dfl)^*^*52v + 2{-l)'^^'^^A*L*BL + 2{-lY^'^^L*A*LB . 

After cancellation this gives 

2J = -2A*L* (lb - {-IY^'^^BL) + 2{-iy^'^^L*A* (lB - {-if^'^^BL 
leading to 

2J =-2(lA-{-IY^'^^AlX(lB-{-IY^'^^Bl\ =-2[L,AYg[L,B], 



\9 ■ 



D 



4.3. A Convergence Theorem. In this section a convergence result will be obtained 
for Lindblad operators with an infinite number of jump operators acting on the algebra 
21(0;). It is worth noticing that the proof of this result is modeled on a similar result found 
in [To]. For each finite set X C £ let Lx be a continuous section of the field 21 such that 

(1) Lx{ijj) G 2tx(i^) and has degree dx, 

(2) the section is covariant, namely rj^^^a {Lx-a{^^"'^)) = Lx{uj), 
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(3) for all natural numbers A^ G N there exists p > such that 

N 

^e^^^sup Y. \\Lx{co)r=:CL<oo, 

m6N n=0 ^'^^ OeX ; diam(X)=m;|X|=n+l 

where diani(X) denotes the diameter of the set X and |X| the number of points 
in X. For convenience, the diameter will be calculate w.r.t. the metric in M"' given 
by |x| = max \xi\] 1 < i < d. Hence, a ball of diameter i? is a cube of side R with 
sides parallel to the canonical basis of M"^. Such a cube has a volume W^. On the 
other hand, only the nonempty X's matter, so that |X| = n + 1 for some natural 
integer n. 

Then P^^^j denotes the operator acting on 2t(a;) defined by 

VxA^) = l{LxiujrLxiuj),A} - i-lf-'^^LxiiorALxiuj) , 
where {A, B} = AB + BA is the anti-commutator. Then the following result holds. 

Theorem 4.4. Let 3" denote the set of finite subsets of C. Then, the operator S^^ = 
^xeJ "^^.^ ^'^ '"'^^^ defined on 'Qlioduj). It is covariant and, if ^ is a continuous section 
of '^loc, its image 2)(0 is a continuous section of 'Ql. In addition its exponential e~*®'^ 
defines a continuous covariant field of Markov semi-groups. 

Proof: 1)- Let A E ^a{uj). Then 'Vx,uiiA) = whenever XflA = 0. Hence the sum over 
X is restricted to those subsets X intersecting A. From the estimate given in eq. (tsl), it 
follows easily that the sum defining 1)i^{A) converges in norm, uniformly w.r.t. u. Then 
the continuity and the covariance are straightforward to check. 
2)- Iterating the definition of T) gives (omitting u) 

T>\A)= Y. -Dx.o.-.oVxM), 

where the family [Xi, ■ ■ ■ ,Xk) of finite subsets of C is submitted to satisfy the compati- 
bility condition Xj fl Aj_i 7^ for 1 < j < fc, whenever Aq = A and Aj = Aj_i U Xj. It is 
easy to check that 

\\Vx{A)\\ <2\\Lxf\\A\\. 

Therefore, the l.h.s. can be estimated by 

k 

\\T)'{A)\\<\\A\\ Y 2'^nil^^.H'- 

By assumption, HLxP < (j^^-pdmm(Xj)^ -^^^ N{mi,ni;--- ]m}^.,nk) denote the number 
of compatible families (Xi, ■ ■ ■ ,Xfc) such that diam(Xj) = nij and \Xj\ = Uj + 1. This 
gives 

\\i:)\A)\\<{2CLf Yl e-P('"i+-+™'=) Y N{mun,;---;mk,nk). 

mi,--- ,nn: ni,--- ,rn. 



24 GEORGE ANDROULAKIS, JEAN BELLISSARD, AND CHRISTIAN SADEL 

In order to estimate N{mi, rii, ■ ■ ■ ; mk,nk), it ought to be remarked that 

(i) l-^l < diam(X)'*, so that Uj < m'j for all j, 

(ii) the maximal number of choices for Xj is obtained by choosing a point in Aj_i then 



dui 



by choosing rij points in a hypercube of side at most rrij. There is at most m- ^ x |Aj_i| 
ways of making this choice. 

(iii) the number of points in Aj_i is at most {|A| +rai + n2 + - ■ ""^j-i} < {|A|+iV(j — 1)} < 
Nj{l + (|A|/A^ — l)/j), by construction. In particular this gives 

N{mi,n,;---;mk,nk)<N''k\ fjmf JJ /"l + ^^1^") . 
Then using (1 + u) < e" and also 1 + 1/2 + ■ ■ ■ + 1/k < 1 + ln(/c) gives 

k 



N{m,, ni;---;mk,nk)<N''k\ el^l/^"^ fcl^l/^-^ J] "^^"^ 






It follows that 



fc=0 ■ k=0 

with K = |A|/A^ - 1 and Ci = 2NCl Em>i m^^e-P'" < cx). Hence, if Ci\t\ < 1, the sum 
defining e~^^{A) converges absolutely and uniformly in a; G f2. Since Ci does not depend 
on the volume A it follows that e~*® is well defined on 21 for t E C such that |t| < Ci^. 



3)- Using the Proposition 4.1, and the previous convergence, it follows that e~*® is CP 
for < t < Ci^. Moreover, since 2)(1) = 0, it is actually CPi. In particular, it is a 
contraction semi-group. Moreover, by construction of the exponential, if s, t G C are such 
that \s\ + |t| < Cf\ then e"'^'*"'"*^® = e~*® e~*®. Therefore if t G IR+ let n be an integer 
such that t/n < Cf ^. Then e~*® = (e~*/"®) is well defined and does not depend upon 
which n has been chosen. And for the same reason it defines a Markov semi-group. 

4)- The continuity and the covariance follow from the definition and the proof will be left 
to the reader. □ 



4.4. Jump Dynamics. The general dissipation operator 2)^^ defined in Section 2.5 is the 



main object of study in this subsection. The main result is Theorem 4.9 



The model defined in Section 2.5 is a specific example of a larger class of models of the 
form 

®-(^) = E (^{^;M^7M> M - i-ir'^'-L;{u)AL,{u)] , (4.4) 

7e3 ^ ^ 

where the following axioms are satisfied. 

Jl- The index set 3, called the set of jumps, is countable and admits a bijective 
involution 7 G H *— ?■ 7 G ^ called time-reversal. Moreover, the translation group Z"^ 
acts on ^ in a bijective way and the action is denoted by 7 i—)- 7 + a by mappings 
commuting with the involution. 
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J2- For each 7 G 3 there is a continuous covariant field of local observables L-y, called 
the jump operators. In particular 

(i) 3 A C £ finite, depending on 7, so that L^{(jj) G ^a{(jj) , Vw G fi, 
(ii) the smallest such A is called the support of 7 and is denoted supp{7} and 
satisfies supp{7 + a} = supp{7} + a, 
(iii) rji^^a {-^7-a(T~"a;)} = L^^u) and L^ , L^-a have the same degree. 

J3- Under the time evolution the jump operators satisfy at{L^) = 6***^^ L^ for all t G M. 
where the e.y : Q ^ M's are continuous functions. 

J4- The jump operators satisfy a (3-KMS condition., namely L* = e~^^''^'^ Ly, in par- 
ticular e^ = —By and Ly has the same degree as L^. 

J5- The following sum converges ^ ^^ r 1 L*(u;)L^(u;) G 21(0;) uniformly with 
respect to uj. 

Remark 4.5. The /3-KMS condition is also called detailed balance in the Physics liter- 
ature (see for instance |12]y'- This axiom is essential to describe the dissipative evolution 
at equilibrium. 

Proposition 4.6. If the jump operators satisfy the assumptions [J1-J5] {except possibly 



the axiom fJ4j), the operator T)^ is well defined by eq. (4-4) on the set of local observables. 
In addition 

of' o 2)^ o a I = 2)(^ , (time-invariance) (4.5) 

Vio,a o '^T—Lo o VZ!a = ®^ ' (space-covariance) (4.6) 

and, if C, E 3^, the field u E Q \-^ 'D^{^{uj)) is continuous. 

Proof: Let A be localized in Aq. Each L^ belongs to the local algebra 2tsupp{7}- 
In particular if supp{7} fl Aq = 0, it follows that A commutes with L*L^ and that 
{-ifi^i'^A^ALy = L*yLyA. Thus, the term |{L;L^, A] - {-lY^-^'^^^AL^ just vanishes. 
Hence, 

7Ga;SUpp{7}nAo7^0 ^ ^ 

Thanks to the axiom [J5], this sum converges in 2t. Thanks to axiom [J2-J3] the time- 
invariance, the covariance condition and the continuity are satisfied. □ 

Proposition 4.7. Let p be the continuous field of (5-KMS states over 21 defined in 



Prop. 3. 1 . If the jump operators satisfy [J1-J5], then p is 'D-invariant and satisfies for all 



finite A C C and all A,Bg ^a{uj) 

p.iA*J)UB))=pU^UArB) = IJ2p- {[L,ico),A];[L,iu),B],) . (4.7) 



2 

7Ga 



Remark 4.8. This shows thatT) acts on the GNS-representation of the ground state as a 
a positive self-adjoint operator. Moreover it can be seen as a generalization of a Laplacian. 
For indeed, the family of (graded) commutators by L^, indexed by •j G 3, can be seen as a 
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gradient and the right hand side of the previous equation looks like a Sobolev norm of the 
type /|VA|2 if A = B. 

Proof: Since there is no confusion, the reference to w G fi will be omitted. Thanks to the 



Leibniz formula (Proposition 4.3), it is sufficient to show that (i) po2) = (2)-invariance 
of p) and (ii) 2) is symmetric. 

1)- Let A G SIa. Since 1D{A) has the same degree as A, it is sufficient to assume that dA = 
0, because p{B) = whenever B has degree one. Then p {1){A)) is the sum of three types 
of terms, namely {l/2)p {L*^L^A), {l/2)p {AL*^L^) and {-!) p {L*^AL^) . Thanks to [J3], 
it follows that L*L^ is invariant by the dynamics at- In particular a_j/3(L*L^) = L*L^. 



Thus, tanks to the /3-KMS condition (eq. (3.10)) 



p(AL;L^) = p (a„,^(L;L^)A) = p {L*^L^A) . 

In particular the ffist two terms are equal. Moreover, the axioms [J3-J4] imply that 
a-rf3{Ly) = e^^-'Ly = e^^^/^Li, and, similarly, c'^^t/^L* = Ly. Hence 

p{L*^AL^) = p{a.,p{L^)L*^A) = p{L*^L^A) . 

Therefore, since the time-reversal 7 H- 7 is a bijection, the sum of these last terms 
compensate the sum of the other terms to give zero. 

2)- To prove that 2) defines a symmetric operator, let A, B be elements in 21a. Without 
loss of generality, it can be assumed that A and B have the same degree, otherwise the 
p-average vanishes and the identity becomes trivial. Then p{A*'D{B)) is a sum of three 
types of terms. The first ones are 

p{A*L*^L^B) = p {{L^L^AyB) . 

The next terms have the form 

p{a*bl;l.,) = p {l;l^a*b) = p {{AL^L^yB) , 

where the middle identity comes from the /3-KMS condition and the invariance of L*L^ 
under the time evolution (axiom [J3]). The last terms are more involved. Using again the 
/3-KMS condition, the identities d^ = ds, d^ = d^ and the axiom [J4], leads to 

(-i)^-^^p(a*l;5l^) = {-if-r'^^e'^'^p {l^a*l;b) = {-if^'^^p (mAL^yB) , 

Since the map 7 i—)- 7 is a bijection the sum of all these terms gives p{A*1){B)) = 
p(2)(A)*fi). D 

Through the GNS-representation, the field of ground-states p defines a continuous trans- 
lation covariant field of Hilbert spaces "Ho; = -^^(2t(u;),paj) [231 US]- ^w is obtained from 
2t(a;) through the inner product {A\B)i_j = pi^{A*B), after taking the quotient by the 
subspace of elements of zero norms and completing. The canonical map from 21(0;) into 
Tiij will be denoted by Cw- Hence 

p^{A*B) = {UA)\UB))., pUA*A) = WUml, ABe 2l(a;). 
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On each element of this field the Proposition 4.7 defines a densely defined field of positive 
quadratic forms as follows 

QM.B) = i5^p4[L,H,A]*[L,(u;),5y . (4.8) 

The following result shows that, as a consequence of the Friedrich extension theorem [57] , 
this form defines a positive self adjoint operator "H^. 

Theorem 4.9. // the jump operators satisfy [J1-J5], the quadratic form Q^j, which is 
densely defined on 7{^, is closahle. Its closure defines a positive self-adjoint operator, 
denoted by T>^, on Ti^. The contraction semi- group e"*^'^ is completely positive. The 
corresponding field of contraction semi-groups is continuous, time-invariant and covariant. 

Proof: (i) Closability. By abuse of notation A will represent here either an element of 
2l(w) or its image in "H^^. Let ||A||^^q denote the norm 

\\A\\Iq = \\A\\1 + Q^{A,A) 

To prove that Q^^ is closable, let {An)n£N be a || • ||^,Q-Cauchy sequence in 'Qliod'-^) such that 
lim„__i.oo ll^nlU = 0. It should be proved that lim„_j.oo ||^n||a;,Q = 0. Since this sequence is 
II ■ ||a;,Q-Cauchy, it follows that Quj{An, An) converges and is therefore uniformly bounded in 
n. In particular, it follows that S{An) = {[L^/{u), A^j^) ^ is Cauchy if seen as an element 
of the Hilbert space H^j ® i'^{3)- Therefore there is ^ = {^j)^^ e U^ ® i'^{3), such that 
S{An) — )■ ^. Hence, given e > there is a finite subset I G 3, such that J^-yi^i ll'C7llw — ^■ 
Now, thanks to [J1-J5], and omitting cu inside L^ 

Puj (l|-^7;^njg| ) = Puj (^^"^7-^7^") + Pw (-^7^n^"-^7J 



_(_i)'^..^A„ (^p^ (a;l;a„l^) + p^ (l;a;l^a„) 



The last two terms of the r.h.s. can be estimated in terms of the first two. Since L^iu) 
is bounded, the first term is bounded by ||L^(a;)||^||y4„||^ which converges to zero. Using 
the axiom [J3-J4], the second term can be written (omitting w), as 

P \I'^A^AnL^j = e '' p (^L^^L^Aj^An) = p (^L^L^A^An) ■ 

Thanks to [J3], it follows that Lt-L^y is invariant by the modular automorphism, so that, 
using the Cauchy-Schwartz inequality, 

p mL^A*^An) < \\AJp {Lt^L^AlAn^L^y^ = ||A„||p ((L^L^)'a;A„)'^' . 

Iterating m-tiines, leads to 

n(T*T A* A \<\\A \\^-'^'""n( (T*T "l^'" A* A 1 < II 4 l|2-2"™||7' W^W A ll^"" 

/^ V 7 7 n "/ — ll^"ll P {{-^J-^lJ ^n^nj ^ ll^nll ll-'^yll ll^njla 

Consequently, as m — )■ 00, 

n ( T* T A* A ^ < II 4 ll^ll 7' l|2 

P \^^J^^^n^n) ^ ll^nll ||-'^7l| ; 
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which is also converging to zero as n — )■ oo. Hence, the finite sum ^^g/P (|[-^7, ^JP) 
vanishes, so that ||^|| < e. Since e can be taken arbitrarily small, it follows that .^ = 0, 
proving that Q^j is closable. 

It follows from the Friedrich extension method that there exists a positive self-adjoint 
operator D^, with domain V = {A e n,,; 3C > 0, \\Quj\\UB,A) < C\\B\\^ , \/B e V.^], 
defined by 

Q^{B,A) = {B\V^A)^, AeV. (4.9) 



(ii) Complete Positivity. Let 3 be finite to begin with. Thanks to the Propositions 4.7 and 



4.1 , it follows that the complete positivity holds if 3 is finite in the algebra 21(0;). Moreover, 



since 2) commutes with the dynamics (Proposition 4.6), the operator T)^ commutes with 
the Modular operator A^^. In particular e~*^" will be completely positive on "H^ (see 



Section 4.1). If S is not finite, then Q^ can be seen as the supremum of a countable 



family of similar forms with 3 finite. Correspondingly there is a non decreasing sequence 



of positive self-adjoint operators I'n,^ converging weakly to D^^. The Lemma 4.10 belows 
shows then that 2?n,aj converges to "D^^ in the strong resolvent sense. This, in turns proves 
that the semi-group e~*^"" converges strongly to e~*^" . In particular, e~*^" is completely 
positive. 

(iii) Covariance and Continuity. The same argument as before shows that the continuity 
of the field a; G f2 H- Q^ implies the strong resolvent continuity of the field a; G f2 H- 2)^^. 
Hence the semi-group e~*^" is also continuous in u. The covariance is a simple consequence 
of Proposition 4.6 □ 



Lemma 4.10. Let A be a positive self-adjoint operator on the Hilbert space %, with 
dense domain V. Let {An)nm be a non decreasing sequence of bounded positive operators 
converging weakly to A onV. Then the resolvent of A^ converges strongly to the resolvent 
of A. 

Proof: Since Q < Ai < A2 <■■■< An <■■■< A, ii follows that the sequence 
of bounded operators (1 + An)~^ converges weakly. Let R be the weak limit. Since 
the sequence is bounded the sequence converges strongly as well. Then the inequahty 
1 > (1 + An)'^ > (1 + A)'^ implies 1 > i? > (1 + A)'^. In particular, 

1 + A>{1 + Afl\l + Any\l + Af'^ > (1 + Afl^Ril + Af/^ > 1 . 

Hence, (1 + Ay/'^R(\ + A)^^"^ is invertible and its inverse satisfies 

(1 + A)-^/^{l + A„)(l + A)-^/^ < (1 + A)-^'^R-\l + A)-^/^ < 1 . 

Since the left hand side of this inequality converges weakly to 1, thus also strongly, it 
follows that R-^ = (1 + A)- Thus s - lim„^oo(l + A^Y^ = (1 + A)'^ . By standard 
arguments, this shows that for 2; G C \ ]R_|_, s — lim„_j.oo(^n — z)~^ = {A — z)~^. □ 

5. The return to equilibrium 
In this section we first justify the introduction of the thermal bath part of the dissipation 



operator in eq. (2.16). Then it will be shown that the dissipation operator £(^, given in 
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eq. (2.17), defines a semigroup on the state space. Furthermore we obtain that the point 



wise hmit of the semigroup at any initial state is equal to the equilibrium state as the 



time parameter tends to infinity (Theorem 5.4). This justifies the title of the present 



Section. Finally, the spectrum of the dissipation operator 2)(^ is examined which allows 
us to conclude that the return to equilibrium as exponentially fast in time. 

If 2) is the generator of a semigroup (e~*®)(>o on a C*-algebra 21, then an element ^ e 21 
is called invariant for the semigroup if A is an eigenvector of e~*® corresponding to the 
eigenvalue 1 for every t > 0. Also, if p is any state on the C*-algebra 21 then p is called 
©-invariant if 

p{e-'^{B)) = p{B) VBG2t,Vt>0. 

The following result states elementary properties of the jump operators which are defined 



by eq. (2.9). It can be obtained by inspection and the proof will be left to the reader. 



Proposition 5.1. The jump operators Lx^y{uj) given by eq. (2.9) & (2.11) have degree 
zero and satisfy the axiom [J1-J5] with 

(i) suppjx-^y} = {x,y}, 

(a) {x — > y} + a = {x + a — )■ y + a} 

(Hi) the time-reversal corresponds to {y -^ x}, 

(%VJ Sx^ty ^ ^y ^x- 

In order to justify the need for the thermal bath part of the dissipation operator, it is 
enough to observe that the kinetic part 1)^"' suffers from an annoying disease: since 
it leaves the number operator invariant, it does not have a unique invariant state. In 



fact, varying the chemical potential /i in (3.13) produces infinitely many invariant states. 



Thus, the thermal bath fixes the chemical potential and the kinetic part of the dissipation 
operator alone cannot drive the electron gas towards equilibrium. This is the content of 
the following result. 



Proposition 5.2. The kinetic part of the dissipation operator, given in eq. (2.13), com- 
mutes with the number operator. In particular its set of invariant state is not reduced to 
a point. 



Proof: The number operator has been defined in eq. (3.4) for finite volume. As for 
the Hamiltonian dynamics it generates an automorphism group which is defined, in the 
infinite volume limit by 

v['^\sxax) = C'^s^ax , z/f'^-'(s^.4) = e'^s^al . 

In particular, Lx-^y{uj) is invariant by this automorphism group. Therefore, like in Propo- 
sition 4.6[ it follows that 

Hence any state generated by a finite volume Hamiltonian of the form 

Fk{uj) -- 
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where P and Q are polynomials, is invariant by S)^!,'". □ 

The following result states elementary properties of the jump operators which are defined 



by eq. (2.14). It can be obtained by inspection and the proof will be left to the reader. 



Proposition 5.3. The jump operators Lx^i,{u) and L^_^a.(u;), given by eq. (2.14) & 



(2.15), have degree one and satisfy the axiom [J1-J5] with 
(i) supp{x —!■*} = {x} = supp{-k — )■ x} 
(a) {x — 7- *} + a = {x + a — )■ *} and similarly for {-k — )■ x} 
(Hi) X ^ -k is time reversed from * — )■ x, 
(iv) e^^^ = e^- fi. 



The following theorem shows the return to equilibrium p^ as defined in Proposition 3.7 

)t>o- 



for the dynamical system (e ^^^ 



Theorem 5.4. The operators 2)* and 2)^ have a unique invariant state given by p^. 
Also, for every state p on the C*-algebra Ql{u), and for every observable A G 2t(ci;), 

limp{e''^A)=p^{A). (5.1) 

t—^oo 



Proof: 1)- Assume that 2)* (A) = 0. Then thanks to Proposition |4.7 it follows that 



= p^iA*^UA)) = E.r.^^Pc. i\K,A]f) + E.r,^:.P- {\[alA]f). Thanks to Propo- 
sition |3.8[ Pi^ is faithful, so that 

[a^, A]g = 0= [4, A]g Vx e C{.s) . 

Then the only elements of 2t(a;) with this property are the multiples of 1. To prove this 
let y4 be a local observable first, namely A G 21a (cu) • Then, using the decomposition into 
monomials A can be written as A = n^B + (1 — nx)B' + a^C + al.C' where 5, -B', C, C G 
2tA\{^}(a;). Thus [a^, A\g = C + aj.{B — B') = implies that C = and B = B'. In much 
the same way [al, A]g = implies C = 0. Thus A = B E 'Qi\\{x}{uj). Using this argument 
inductively on every point of A shows that A is a multiple of 1. 

If now A G 2t(a;), then for all We > there is a finite set A^ and A^ G 2tA^(a;) such 
that \\A — Ae\\ < e/3. For any finite subset A of £ let Ea : 2l(a;) — )• 21a (cu) be the 
canonical projection. The commutation rule above implies that Ea(v4) also commutes 
with all the ax,al. for x G A, showing that there is c(A) G C such that Ea(v4) = c(A)l. 
ForAD A^, ||Ea(A-A,)|| = ||c(A)l-A,|| < e/3. It follows that (i) |c(A)| < \\A\\ and (ii) 
|c(A)— c(A')| < 26/3 for A, A' D A^. Thus it is a Cauchy sequence which converges to c G C 
as the volume tends to infinity. Hence ||y4 — cl|| < ||y4 — A^H + lc — c(A)| + ||y4e — c(A)l|| < e. 
This shows that A = cl. 
2)- Assume that D(A) = 0. Then 

= pUA*^UA)) = p^{A*^'J''iA)) + pUA*KiA)) . 
By Proposition Unapplied to 2)f;^" we obtain that p^(A*D^*"(A)) > 0. Thus 

= pUA*^t'''m=pUA*^UA)). 
Hence by 1), A is a multiple of the identity. 
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3)- By construction, both S),^ and 2)* leave each 21a (w) invariant. Therefore e~*®" and 
e~*®'^ are well defined on each of the 2tA(w)'s and define Markov semi-groups (Proposi- 
tion 4.1). Since the multiples of the identity are the only invariant observables of the 
semigroups, 1 is a simple eigenvalue of these semi-groups. Consequently limi_>.oo e~*^"(y4) 
and limi_5.oo e~*®"(y4) exist for each A G Qiiodoo) and, (since these limits are invariant ob- 
servable for the corresponding semigroups), these limits are multiples m^(y4) and m*(y4) 
respectively of the identity. Since the semi-groups are contractions the result extends by 
continuit y to all elements of 21(0;) 



Theorem 4.4 allows us to define e *®'^. For every A G 21^, this leads to 



lime-*®"(A) =m^(A)l and lim e~*®"(A) = m^(A)l . (5.2) 

t—^oo t— >oo 



4)- By eq. (5.2) it follows that m^ and m* are states on 2l(a;). If now p is a S^-invariant 



state and A gQI^, then 

mUA)= limp(e-*^"(A)) = p{A) . 

Similarly, if p is a S)*-invariant state then m* = p. 



5)- Thanks to Proposition 4.7, it follows that p^ is both 2)* and ©-invariant, thus p^^ = m^ 



and puj = m*. Therefore p^j is the unique invariant state for "D^^ and 2)* 



6)- Now let p be any state on 01^. By eq. (5.2), 

lim p (e-*^-(A)) = p(m^(A)l) = m,(A) = p^A) 

which shows the return to equilibrium p^j. □ 



As explained in Section [4^ the dissipation operators, defined previously, define self-adjoint 
positive operators on the Hilbert space of the GNS representation. Some spectral prop- 
erties are given by the following result 

Theorem 5.5. (i) IfD"^ and D^^ denote the corresponding positive self-adjoint operators 
acting on the GNS representation, then both admit as a simple eigenvalue and both have 
a positive gap separating zero from the rest of the spectrum bounded from below by r^/2. 

(a) Let IC^ be the closed subspace of Ti^ generated by vectors in vr^(lt^)^a; which are 
orthogonal to ^^ . Then both 2)* and !D^ leave /C^^ invariant and their restriction to /C,^ is 
bounded below by T^. 
(Hi) The operator 2)* has a pure point spectrum. 

Remark 5.6. This result shows that the Markov semi-group generated by 2)^ converges 
to equilibrium exponentially fast with a lifetime given by the inverse o/F^. 

Proof: (i) In the Hilbert space Tiuj of the GNS representation, monomials in the 
annihilation-creation operators make up a total set. This helps creating an orthonor- 
mal basis in this Hilbert space. First, for x G £ let b^i^oj), &J,(w) be defined by 

b^iu) = e-'^('^-^)^/2 s, a, , bUu) = e-^^^-'^)^/^ s, al . (5.3) 
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It is easy to check that both elements define unit vectors in "H^, and that they are 
orthogonal to each other. For X = (a;i,X2 ■ ■ ■ iXm) G C{(jo)^^ let hx{oo),h*-^{uj) G l-iu) be 
defined by (omitting the reference to u) 

bx = hA,---b.^, b^x = bljl^_,---bl. (5.4) 

If two components of X are equal, then bx = 0. Thus the X's will be restricted to the set 
of elements in C{u)^"' made of distinct points. If the order in presenting the points of X 
is changed, then bx changes sign according to the signature of the permutation. Hence, 
modulo a sign, bx depends only upon the set {a:i, 0:2, ■ ■ ■ , Xm} C C{u). In much the same 
way the vectors crx{oo) will be defined as follows 

a,{uj) = s, (e''/2(^-^)n, - e-/^/2(^-M)(i _ ^j) ^ ^^ = ^^^^^^ . . . ^^^ (5.5) 

Here, the ordering of points is irrelevant since the n^^s commute. Before continuing the 
proof we need the following lemma and corollary. 

Lemma 5.7. The vectors in "H^ given by CuiibxbyO'z), where X, Y, Z vary among the set 
of triplets of three disjoint finite subsets of C{u), including the empty set, make up an 
orthonormal basis ofHi^. 

Proof: First, these elements define all possible monomials in 2t(a;) up to a scalar 
multiplication. Consequently they generate a dense subspace. If it is proved that these 
vectors make up an orthonormal family, then they make up a Hilbert basis. Second, due 
to the factorization property of the equihbrium state p^^ and since X, Y, Z are disjoint, it 
follows that 

Puj ({bxbYO-zTb^bYO-z) = Puj{bxb^x)P^{b\.bY)pu' '''^' 



Again the factorization property of p^^ and the commutation rules gives 

n 

P.{b'^Yby) = X{e-^\'^^-^\'')p^{ny^) = l. 

i=l 

A similar result applies for Pujibxbx)- K also applies to P(^(o-|) once it is remarked that 
'^x — ^x [e^'^'^'^'^^^rix + e~^*^^"=~^)(l — Ux))- Hence each of these vectors is normalized. 

Let now (X, Y, Z) be three disjoint finite subsets of £, and let (X', Y\ Z') be three other 
ones, such that these triples are distinct. Then without loss of generality, it can be assumed 
that X U F U Z 7^ 0. Assume first that there is a point x e {X\JY \J Z)\{X' \JY' \J Z'). 

If a; G X, then the inner product p^; ( {bxbvo'z)* bxibyO'z' ) will have a factor Puj{bx) = 0, 

if a; G y, the factor will be p^ip'^x) — 0) whereas if a; G 2' it will be p^icrx) = 0. Hence the 
two vectors are orthogonal. li X UY U Z = X' UY' U Z', then either X ^ X',Y ^Y' 
OT Z ^ Z'. In the first case there is a; G X \ X'. Then x & Y' U Z'. If a; G Y', then the 
inner product get a factor Puoipxbx) = 0, whereas if a; G Z' the factor is Puj{bxCrx), which 
vanishes as well. The same argument holds in all other cases leading to the result. □ 
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Corollary 5.8. The family {Cca{<^z) ', ^ ^ Z G C{uj) , Z finite} is an orthonormal basis 

Proof: Since elements of IC^j are orthogonal to ^^j = Ca;(l) they are generated by linear 
combination of the bxbyO'z^s with XUYUZ ^ i/}. Since they come from the commutative 
sub-algebra C^j they are generated by the operators Sxri^, namely by the ax's only. □ 



Proof of Theorem 5.5 (continued): Let 'Dx,i, be the operator on 2t(a;) defined by 

1)xA^) = 2 (4«x-4 + Aala^) - 4.cr(A)a^ . (5.6) 

In much the same way, let Tt^^^x be defined by 

'D^,x{A) = - {axalA + /la^4) - a^a(A)4 . (5.7) 

It follows that a A = AxB with Ax e 2t{^} and B ^ 'Ql{x} then 1)xA^) = '^xA^x)B and 
the same holds for 2)^^x- Since Ax can only be a linear combination of 1 , bx , b^ , 0"^;, it is 
enough to consider each of these cases. This leads to 

2 5 '^x,-ky'Jx) 2 



2),,,(1) = 0, ^xAbx) = ^, ^x,M) = f, 2),,,(a,) = 2cosh(f(e,-/i))n, 



D,,,.(l) = 0, 'D.^bx) = |, D.,.(6D = |, ^.,xi^x) = -2 cosh(f (e, - /i))(l - nj . 

Since 2)* = X]xe£(s) ^x^-k'^x,* + r*^.a;2)*,x it follows immediately that bx and 6j_. are both 
eigenvectors of !l)* for the common eigenvalue 

7. = %^ + %^ = y(l + e-'^l--l) . (5.8) 

This means Cuj{bx) and C^l^D are eigenvectors of 2)* with eigenvalue ■jx- Similarly, an 
elementary calculation shows that 

^UUb^xbyaz)) = lxy,z Ub^xbr^rz) , (5.9) 

where 

7x,Y,z= Yl 7. + 25^7. • (5.10) 

xexuY zez 

Hence, D* has pure point spectrum. In addition, 7x,y,z vanishes if and only if X = y = 
Z = 0. This shows that C,ui = Caj(l) is the only eigenvector with eigenvalues 0. It shows 
also that /C^^ is invariant by D* and, thanks to the Corollary 5.8 , the restriction of D* on 
]Ci^ is bounded below by T^ (uniformly in /3,/i and u). 

Let ^x,y{A) = \{Ll^yLx^yA + ALl^yLx^y) - Ll^yALx^y, then one has 1)x,y{nx) = 
Tx^yUxil-ny), ^x,y{ny) = -Tx^yUxil-ny) and ©^^^^(n^) = for 2; ^ {x,y}. Hence, T)^ 
preserves the sub-algebra ft^j. Moreover 2)a;(l) = 0, showing that D^; leaves the subspace 
K-oj invariant. As D^; > D^ it follows that is a simple eigenvalue of T)^ and that Ti^ has 
at least the same gap as D*. □ 
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Remark 5.9. // F^ is not constant but depends on e^ and /3 as specified by Remark 2^ 



then (5.9) and (5.10) are still correct. Only the values of'jx would change correspondingly 



as one has to replace T^ by T^{ex,P) in (5.8). 



By the assumptions of Remark 2.^ the infimum of this constant for all e^ G A and all 



13 G (0, oo) is positive, and therefore D* and T>^ will still have a spectral gap (uniformly 



in j3, /i and u) as asserted in Theorem 5.5 
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